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Abstract 

In this article we compare the characteristics of two types of the intermit- 
tent behavior (type-I intermittency in the presence of noise and eyelet in- 
termittency taking place in the vicinity of the chaotic phase synchronization 
boundary) supposed hitherto to be different phenomena. We show that these 
effects are the same type of dynamics observed under different conditions. 
The correctness of our conclusion is confirmed by the consideration of dif- 
ferent sample systems, such as quadratic map, Van der Pol oscillator and 
Rossler system. Consideration of the problem concerning the upper bound- 
ary of the intermittent behavior also confirms the validity of the statement 
on the equivalence of type-I intermittency in the presence of noise and eyelet 
intermittency observed in the onset of phase synchronization. 
Keywords: fluctuation phenomena, random processes, noise, 
synchronization, chaotic oscillators, dynamical system, intermittency 
PACS: 05.45.Xt, 05.45.Tp, 05.40.-a 



* Corresponding Author 
Email address: moskalenkoOnonlin. sgu.ru (Olga I. Moskalenko) 



Preprint submitted to Physics Letters A 



February 19, 2013 



Introduction 



Intermittency is well-known to be an ubiquitous phenomenon in nonlin- 
ear science. Its arousal and main statistical properties have been studied 
and characterized already since long time ago, and different types of in- 
termittency have been classified as types I-III intermittencies l|, |2| , on-off 
intermittency js, 4|, eyelet intermittency js-?] and ring intermittency jsj. 

Despite of some similarity (the presence of two different regimes alternat- 
ing suddenly with each other in the time series), every type of intermittency 
is governed by its own certain mechanisms and characteristics of the inter- 
mittent behavior (such as the dependence of the mean length of the lami- 
nar phases on the control parameter, the distribution of the laminar phase 
lengths, etc.) of different intermittency types are distinct. There are no 
doubts that different types of intermittent behavior may take place in a wide 
spectrum of systems, including cases of practical interest for applications in 
radio engineering, medical, physiological, and other applied sciences. 

This article is devoted to the comparison of characteristics of type-I in- 
termittency in the presence of noise and eyelet intermittency taking place in 
the vicinity of the phase synchronization boundary. Although these types of 
intermittency are known to be characterized by different theoretical laws, we 
show here for the first time that these two types of the intermittent behavior 
considered hitherto as different phenomena are, in fact, the same type of the 
system dynamics. 

The structure of the paper is the following. In Sec. [1] we give the brief 
theoretical data concerning both the type-I intermittency with noise and 
eyelet intermittency observed in the vicinity of the phase synchronization 
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boundary as well as arguments confirming the equivalence of these types 
of the dynamics. The next Sections |2H3] aim to verify the statement given 
in the Sec. [1] by means of numerical simulations of the dynamics of several 
model systems such as a quadratic map, Rossler oscillators, etc. Eventually, 
in Sec. H] we discuss the problem of the upper boundary of the intermittent 
behavior. The final conclusions are given in Sec. [51 

1. Relation between type-I intermittency with noise and eyelet in- 
termittency 

First, we consider briefly both eyelet intermittency in the vicinity of the 
phase synchronization boundary and type-I intermittency in the presence of 
noise following conceptions accepted generally. The main arguments con- 
firming equivalence of these types of the intermittent behavior are given 
afterwards. 

1.1. Type-I intermittency with noise 

The intermittent behavior of type-I is known to be observed below the 
saddle-node bifurcation point, with the mean length of laminar phases T 
being inversely proportional to the square root of the criticality parameter 
{Sc - e), i.e. 

T~(e,-^)-V^ (1) 

where e is the control parameter and Ec is its bifurcation value corresponding 
to the bifurcation point j^. The influence of noise on the system results 



in the transformation of characteristics of intermittency 
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121 1, with the 



intermittent behavior being observed in this case both below and above the 
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saddle-node bifurcation point Sc- In the supercritical region 12| of the control 
parameter values (i.e., above the point of bifurcation, e > Ec) the mean length 
T of the laminar phases is given by 



1 



; exp 



(2) 



ky/e — Ec " V 3D 
where k = const, D is the intensity of a delta-correlated white noise ^„ 
[(^n) = 0, (^n^m) = D6{n — m)], with Equation ([2]) being applicable in the 
region 

^2/3 < |£ _ < 1 (3) 



of the control parameter plane 10|,ll3|. In this region the criticality parameter 



{e — Ec) is large enough and, therefore, the approximate equation 



InT = B(e - Er?^"^ - Ink 



(4) 



(where B = const) is used typically (see llj for detail) instead of (|2]). In 
turn, the distribution p(r) of the laminar phase lengths r is governed by the 
exponential law 12 1 

p{T)=T-'exp{-T/T). (5) 

1.2. Eyelet intermittency 

For the chaotic systems in the vicinity of the phase synchronization 
boundary (if the natural frequencies of oscillator and external signal are de- 
tuned slightly) two types of the intermittent behavior and, correspondingly, 
two critical values are reported to exist {sl, O, [l^ . Below the boundary of the 
phase synchronization regime the dynamics of the phase difference A(/?(t) 
features time intervals of the phase synchronized motion (laminar phases) 
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persistently and intermittently interrupted by sudden phase slips (turbulent 
phases) during which the value of |A(y9(t)| jumps up by 27r. For two coupled 
chaotic systems there are two values of the coupling strength ei < £2 be- 
ing the characteristic points which are considered to separate the different 
types of the dynamics. Below the coupling strength value ei the type-I in- 
termittency is observed, with the power law T ~ (ei — e:)^^/^ taking place 
for the mean length of the laminar phases, whereas above the critical point 
£2 the phase synchronization regime is revealed. For the coupling strength 
e G {61,62) the super-long laminar behavior (the so called "eyelet intermit- 
tency") should be detected. For eyelet intermittency (see, e.g. js], O]) the 
dependence of the mean length T of the laminar phases on the criticality 
parameter is expected to follow the law 

r~exp(K(£2 (6) 

or 

Hl/T)=Co-C,{62-6)-'/', (7) 



(co, ci and k are the constants) given for the first time in 15|] for the transient 
statistics near the unstable-unstable pair bifurcation point. The analytical 
form of the distribution of the laminar phase lengths has not been reported 
anywhere hitherto for eyelet intermittency. 

The theoretical explanation of the eyelet intermittency phenomenon is 
based on the boundary crisis of the synchronous attractor caused by the 
unstable-unstable bifurcation when the saddle periodic orbit and repeller pe- 
riodic orbit join and disappear . This type of the intermittent behavior 



has been observed both in the numerical calculations 



HE 



6| and experimental 



studies [7[ for the different nonlinear systems, including Rossler oscillators. 



1.3. Theory of equivalence of the considered types of behavior 

Although type-I intermittency with noise and eyelet intermittency taking 
place in the vicinity of the chaotic phase synchronization onset seem to be 
different phenomena, they are really the same type of the dynamics observed 
under different conditions. The difference between these types of the inter- 
mittent behavior is only in the character of the external signal. In case of 
the type-I intermittency with noise the stochastic signal influences on the 
system, while in the case of eyelet intermittency the signal of chaotic dy- 
namical system is used to drive the response chaotic oscillator. At the same 
time, the core mechanism governed the system behavior (the motion in the 
vicinity of the bifurcation point disturbed by the stochastic or deterministic 
perturbations) is the same in both cases. To emphasize the weak difference 
in the character of the driving signal we shall further use the terms "type-I 
intermittency with noise" and "eyelet intermittency" despite of the fact of 
the equivalence of these types of the intermittent behavior. 

Indeed, the phenomena observed near the synchronization boundary for 
periodic systems whose motion is perturbed by noise (in other words, the 
behavior in the vicinity of the saddle-node bifurcation perturbed by noise) 
have been shown recently to be the same as for chaotic oscillators in the 
vicinity of the phase synchronization boundary 
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16|. Thus, both 



for two coupled chaotic Rossler systems and driven Van der Pol oscillator the 



same scenarios of the synchronous regime destruction have been revealed |16 |. 
Moreover, for two coupled Rossler systems the behavior of the conditional 
Lyapunov exponent in the vicinity of the onset of the phase synchroniza- 
tion regime is governed by the same laws as in the case of the driven Van 
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der Pol oscillator in the presence of noise 13||. Additionally, when the tur- 



bulent phase begins the phase trajectory demonstrates motion being close 
to periodic both for the eyelet intermittency observed in the vicinity of the 
phase synchronization boundary (see [l^) and for type-I intermittency with 
noise. Finally, the repeller and saddle periodic orbits of the same period in 
the vicinity of the parameter region correspondin g to the intermittent be- 
havior tend to coalesce with each other (see, e.g. [slJl^) for both these types 
of the intermittent behavior. Obviously, if the phenomena observed near 
the saddle-node bifurcation point for the systems whose motion is perturbed 
by noise are the same as for chaotic oscillators in the vicinity of the phase 
synchronization onset, one can expect that the intermittent behavior of two 
coupled chaotic oscillators near the phase synchronization boundary (eye- 
let intermittency) is also exactly the same as intermittency of type-I in the 
presence of noise in the supercritical region. 

So, if type-I intermittency with noise and eyelet intermittency taking 
place in the vicinity of the chaotic phase synchronization onset are the same 
type of the system dynamics, the theoretical equations (jlj) and ([7]) obtained 
for these types of the intermittent behavior are the approximate expressions 
being the different forms of Eq. describing the dependence of the mean 
length of the laminar phases on the criticality parameter. Therefore, Eq. ([7]) 
can be deduced from Eq. (j4]) and vice versa. As a consequence, the coeffi- 
cients B, k and cq, Ci in (jlj) and ([7]) are related with each other. Obviously, 
the mean length of the laminar phases must obey Eq. (jl]) and Eq. ([7]) simul- 
taneously, independently whether the system behavior is classified as eyelet 
intermittency or type-I intermittency with noise. Additionally, the laminar 
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phase length distribution for the considered type of behavior must satisfy the 
exponential law ([5]). 

The intermittent behavior under study is considered in the coupling 
strength range Ec < s < 62- In the case of the system driven by external 
noise (type-I intermittency with noise) the lower boundary value Sc corre- 
sponds to the saddle-node bifurcation point when external noise is switched 
off. For the dynamical systems demonstrating the chaotic behavior (eyelet 
intermittency) the lower boundary Sc may be found, e.g., in the way described 



m 



13| . As far as the choice of the upper boundary value £2 is concerned. 



this subject is discussed in detail in Sec. HJof this paper both for chaotic and 
stochastic external signals. 

To find the relationship between coefficients in and ([7]) we introduce 
the auxiliary variable ^ = {s-z — e)"^^^ and expand ln(l/T) determined by 
Eq. (j4]) (type-I intermittency with noise) into Taylor series in the vicinity of 
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the point = {^2 — ^o) ^^^j where Ec < £o < ^2, i-e. 
ln(l/r) = {\nk-B{eo-e,f/')- 



m /c 



(8) 



2-\/^o ~ 



35 (£2 - e^o)^ (Seo - £^2 - 4£c) , ^_l/2 

l£^2-£^J 



Having neglected the term o((^ — Co)^) in dH]) one can write Eq. (|8]) in the 
form 

ln(l/T) = Co - ci(£2 - + C2(£2 - ey\ (9) 

Having required 

35 (£2 - £0)^ (4^0 - £2 - 3£c) „ „s 

= 2VS^. 



we obtain, that 



£2 + 3£c ^ ^ /..v. 

eo = ^ , ec<eo< 82 (H) 



and, therefore, equation ([9]) describing the dependence of the mean length 
T of the laminar phases on the criticality parameter for type-I intermittency 
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ln(l//i) 



-10 



0.48 



I ' I ' I ' I 
0.52 0.56 0.60 (s2-s)-i/2 



Figure 1: The dependence /i(e) = k^^ cxp(Be'^^^) simulating the theoretical law (|4]) and 
its approximation by the curve f2{x) = co — ci(e — £2)^^^^ corresponding to law ([7]) for 
eyelet intermittency. The parameter values are = 23, B = 3.33, cq = 23.5, ci = 51.9, 
X2 — 4.0. The points corresponding to fi{x) are shown by symbols •. The theoretical 
law ([7]) for eyelet intermittency is shown by the solid line 

with noise coincides exactly with Eq. ([7]) corresponding to eyelet intermit- 
tency. Correspondingly, in terms of Eq. ( ITTl) the relationship between coeffi- 
cients B, k and cq, ci in (jl]) and ([7]) is the following 

Co =\nk + B{e2-e,Y/^ 

(12) 

Fig. [1] illustrates the relationship of two theoretical laws (jlj) and (JTj) in the 
region e > Ec- Here function fi{x) = exp{Be^^'^) simulates the theoretical 
law (jl]) (the critical point is supposed to be Ec = 0), whereas the curve 
f2ix) = Cq — 01(6—62)"^^'^ corresponds to law ([7]) for eyelet intermittency. The 
value of coefficients B, k, cq and Ci have been selected according to Eq. f[T^ . 
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One can see that in the region of the study both curves coincide with each 
other. It means that the mean length of the laminar phases obeys Eq. (jl]) 
and Eq. ([7]) simultaneously, independently whether the system behavior is 
classified as eyelet intermittency or type-I intermittency with noise. 

2. Numerical verifications 

To confirm the concept of the equivalence of intermittencies being the 
subject of this study we consider several examples of the intermittent be- 
havior classified both as eyelet intermittency taking place in the vicinity of 
the phase synchronization onset (two coupled Rossler systems) and type-I 
intermittency with noise (quadratic map and driven Van der Pol oscillator). 

2.1. Two coupled Rossler systems 

As we have mentioned above, the intermittent behavior of two coupled 
chaotic oscillators in the vicinity of the phase synchronization boundary is 



classified traditionally as eyelet intermittency (5|, [g 



Nevertheless, the 

behavior of two coupled Rossler oscillators close to the phase synchronization 
onset was considered from the point of view of type-I intermittency with 
noise for the first time in ll|, whereas the same dynamics from the position 
of eyelet intermittency was studied in ^. According to different works the 



mean length of laminar phases happens to satisfy both Eq. (j4]) (Ref. ll| ) 



and Eq. ([7]) (Ref. [6|). Recently [17| the distribution of the laminar phase 
lengths has been found to obey the exponential law corresponding to 
type-I intermittency with noise. To give the complete picture we replicate 
the consideration of two coupled Rossler systems near the onset of the phase 
synchronization regime for the different type of coupling between oscillators 
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and another set of the control parameter values and show that the observed 
intermittent behavior may be classified both as eyelet intermittency and type- 
I intermittency with noise. 

The system under study is represented by a pair of unidirectionally cou- 
pled Rossler systems, whose equations read as 

id = -^dVd - Zd, 
Vd = ^^dXd + aVd, 

Zd=P + Zd{Xd - c), 

(13) 

ijr = OOrXr + ttl/r, 

Zy I Zy (yOC y ^ 

where {xd,yd,Zd) [{xr,yr, Zr)] are the Cartesian coordinates of the drive [the 
response] oscillator, dots stand for temporal derivatives, and a is a parameter 
ruling the coupling strength. The other control parameters of Eq. fll3p have 
been set to a = 0.15, p = 0.2, c = 10.0, in analogy with our previous 



studies 
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19| . The cj^-parameter (representing the natural frequency of the 
response system) has been selected to be Ur = 0.95; the analogous parameter 
for the drive system has been fixed to Ud = 0.93. For such a choice of the 
control parameter values, both chaotic attractors of the drive and response 
systems are phase coherent. The instantaneous phase of the chaotic signals 
(f{t) can be therefore introduced in the traditional way as the rotation angle 
(pd,r = s^T^ctan{yd,r/xd,r) on the projection plane {x,y) of each system. 

In Fig. [2] one and the same result of the numerical simulation of two cou- 
pled Rossler systems f[T^ is shown in different ways to compare obtained 
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Figure 2: The points obtained numerically for two unidirectionally coupled Rossler os- 
cillators ([T3|) are shown by symbols The theoretical laws ([7]) and (j4]) are shown by 
the solid lines, (a) Eyelet intermittency: the dependence of ln(l/T) on the parameter 
((72 — f)^^^^; = 0.042, Co = —0.75, ci = 0.55. (b) Type-I intermittency with noise: 
the dependence of the mean laminar phase length T on the parameter {a — ac)^^^ , with 
the ordinate axis being shown in the logarithmic scale; CTc = 0.0345, B = 1.01 x 10^, 
k = 6.68 X 10"'' 

data with the analytical predictions (JTj) and (jlj) for eyelet intermittency tak- 
ing place near the phase synchronization boundary (Fig. [2], a) and type-I 
intermittency with noise (Fig. |2l b), respectively. The dependence of T on 
((72 — c) is shown in the whole range of the coupling parameter strength val- 
ues (Fig. [2], a) to make evident the deviation of numerically obtained data 
from law ([7j) far away from the onset of the phase synchronization. The cou- 
pling strength a plays the role of the control parameter. The critical point 
(T2 ~ 0.042 relates to the onset of the phase synchronization regime in two 
coupled Rossler systems. The point cXc ~ 0.0345 used in ([2]) and (jl]) corre- 
sponds to the saddle-node bifurcation point if the chaotic dynamics being 
the analog of noise could be switched off. The value of this point has been 
found from the dependence of the zero conditional Lyapunov exponent on 
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the coupling strength (see for detail jl3|). 



One can see, that the intermittent behavior of two coupled Rossler sys- 
tems may be treated both as eyelet and noised type-I intermittency with the 
excellent agreement between numerical data and theoretical curve in both 
cases. Moreover, the coefficients cq, ci and B, k of the theoretical equa- 
tions ([7]) and (jl]) agree very well with each other according to Eq. (fT2|) . It 
allows us to state that both these effects are the same type of the system 
dynamics. Nevertheless, to be totaly convinced of the correctness of our 
decision we have to consider other examples of the intermittent behavior 
classified traditionally (contrary to the previous case of two coupled Rossler 
systems) as type-I intermittency with noise. 

2.2. Driven Van der Pol oscillator with noise 

The second sample dynamical system to be considered is Van der Pol 
oscillator 

X - {\- x'^)x + X = Asm{uet) + D^{t) (14) 

driven by the external harmonic signal with the amplitude A and frequency 
Ue with the added stochastic term DC,(t). The values of the control param- 
eters have been selected as A = 0.1, Ue = 0.98. For the selected values of 
the control parameters and D = the dynamics of the driven Van der Pol 
oscillator becomes synchronized when A = Ac = 0.0238 that corresponds to 
the saddle-node bifurcation. The probability density of the random variable 

m is 

MO = ^expf-i^-V (15) 



27ra V 2a2 
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Figure 3: The points obtained numerically for driven Van der Pol oscillator with the 
stochastic force (|14p are shown by symbols The theoretical laws (0 and ^ are 

shown by the solid lines, (a) Eyelet intermittency: the dependence of ln(l/r) on the 
parameter {A2 — A)^^/^; A2 = 0.029, cq = 1.24, ci = 0.62. (b) Type-I intermittency with 
noise: the dependence of the mean laminar phase length T on the parameter (A — Ac)^^'^\ 
Ac = 0.0238, B = 2.37 x lO"*, k = 4.76 x IQ-* 

where = 1. To integrate Eq. (fT4|) the one-step Euler method has been 
used with time step h = 5 x 10~^, the value of the noise intensity has been 
fixed a.s D = 1. 

On the one hand, as it has been discussed above, the intermittent be- 
havior in this case have to be classified as type-I intermittency with noise. 
The corresponding dependence of the mean length of laminar phases on 
the criticality parameter {A — Ac) is shown in Fig. [3l b. If the amplitude 
A of the external signal exceeds the critical value A^ the exponential law 
T ~ exp(a(y4 — Ac)^^"^) is expected to be observed. To make this law evident 
the abscissa in Fig. [3l b has been selected in the {A — ylc)^/^-scale and the 
ordinate axis T is shown in the logarithmic scale. One can see again the 
excellent agreement between the numerically calculated data and theoretical 
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prediction d?}. The distribution of the lengths of the laminar phases p(t ) 



obtained for A > Ac also confirms the theoretical curve ([5]), see Fig. 7 in [12|. 

On the other hand, trying to choose the corresponding values of A2 for 
the driven Van der Pol oscillator f lT4|) one can find out that the intermittent 
behavior of this system also may be identified as eyelet intermittency. Indeed, 
in Fig. 131 a one can see a very good agreement between the numerically 
obtained mean length T of the laminar phases for the different values of 
the coupling parameter and theoretical law ([7]) corresponding to the eyelet 
intermittency. Note also, that for the well chosen values of A2 the dependence 
T{A2 — A) in the axes ((^2 - A)-^/^ln(l/T)) behaves in the same way as 
the corresponding function T{a2 — cr) in the axes (((T2 — o")^^^^; for 
two coupled Rossler systems f lTS]) . Again, as well as for two coupled Rossler 
oscillators, the coefficients cq, ci and B, k of the theoretical equations ([7]) 
and (jl]) agree very well with each other according to Eq. (fT2|) . 

2.3. Quadratic map with stochastic force 
The next example is the quadratic map 

Xn+i = x"^ + X - e + D$,n, modi, (16) 

where the operation of "mod 1" is used to provide the return of the sys- 
tem in the vicinity of the point x = 0, A = 0.25 and the probability density 
of the stochastic variable ^ is distributed uniformly throughout the interval 
^ G [—1, 1]. If the intensity of noise D is equal to zero the saddle-node bifur- 
cation is observed for e = 0. The intermittent behavior of type-I is observed 
for e < 0, whereas the stable fixed point takes place for e > 0. Having added 
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Figure 4: The points obtained numerically for quadratic map (jl6p are shown by sym- 
bols The theoretical laws ([7]) and (jH) are shown by the solid lines, (a) Eyelet 
intermittency: the dependence of ln(l/r) on the parameter [e2 — s)^^^"^', £2 = 10^^, 
Co = 7.9, ci — 6.45 X 10~^. (6) Type-I intermittency with noise: the dependence of the 
mean laminar phase length T on the parameter e'^/^; = 0, D — 10^^, B = 6.62 x 10^'*, 
k = 2.35 X 10-6 

the stochastic force {D > 0) in ( |T6l) we suppose that the intermittent behav- 
ior must be also observed in the area of the positive values of the criticality 
parameter e, with the mean length of the laminar phases obeying law (j3]). 

Although in this case we deal with type-I intermittency with noise, the 
numerically obtained points corresponding to the mean length of laminar 
phases are approximated successfully both by Eq. ([7]) and (j4]) (see Fig. H]), 
with the coefficients cq, ci and B, k of the theoretical equations (I7j) and (jlj) 
agreeing with each other according to Eq. (fT2l) . These findings confirm our 
statement about identity of the considered types of the intermittent behavior. 

So, having studied the intermittent behavior of different systems which 
(based on the prior knowledge) should be classified either eyelet intermittency 
in the vicinity of the phase synchronization boundary or type-I intermittency 
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with the presence of noise, we can conclude that the obtained characteristics 
are exactly the same in all cases described above. Two next sections are 
devoted to the consideration of another systems to give the additional proofs 
of the correctness of the introduced concept. 

3. Van der Pol oscillator driven by the chaotic signal 

In this section we consider Van der Pol oscillator driven by the chaotic 
signal of Rossler system 

id = a{-ujy - z), 

yd = a{ujx + ay), ^^^^ 

Zd = a{p + z{x - c)), 

ii — (A — u'^)u + u = e{Dy — u), 

where a = 0.15, p = 0.2, c = 10, A = 0.1, u = 0.9689 are the control pa- 
rameters. The auxiliary parameters a = 0.99 and D = 0.0664 alter the 
characteristics (the amplitude and main frequency) of the chaotic signal in- 
fluencing on Van der Pol oscillator. 

From the formal point of view the behavior of system f|T7|) can be classi- 
fied neither eyelet intermittency nor type-I intermittency with noise. Indeed, 
since the response oscillator is periodic there are no unstable periodic orbits 
embedded into its attractor to be synchronized, therefore, the system dynam- 
ics can not be considered as eyelet intermittency. Alternatively, due to the 
presence of chaotic perturbations there is no pure saddle-node bifurcation in 
this system to say about type-I intermittency. Nevertheless, it is intuitively 
clear that this example is nearly related to all cases considered above and 
one can expect to observe here the same type of intermittency as before. 
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Figure 5: The points obtained numerically for Van der Pol oscillator driven by chaotic 
signal of Rossler system ([TT]) are shown by symbols "•" . The theoretical laws ([7]) and 
(|4l) arc shown by the solid lines, (a) Eyelet intermittency: the dependence of ln(l/T) 
on the parameter (£3 - e)-^/^. = 0.023, cq = -3.50, ci 2.46 x 10"^ (b) Type-I 
intermittency with noise: the dependence of the mean laminar phase length T on the 
parameter (e - 6^)^/^ Ec = 0.0185, B ^ 1.25 x 10"^, k 6.90 x lO"'^ 

Fig. [5] makes this statement evident. Indeed, the numerically obtained 
data obey both laws (JTj) and (jlj), with the coefficients cq, ci and 5, /c of the 
theoretical equations ([7]) and (jl]) agreeing with each other in accordance with 
Eq. f lT2|) . Additionally, the distribution of the lengths of the laminar phases 



follows the exponential law 



20|, that allows us to say that we deal here with 



the same type of the dynamics as in the cases of quadratic map f ITB]) . driven 
Van der Pol oscillator (fT4|) and two coupled Rossler systems ( !T3|l considered 
above. 



4. Upper boundary of the intermittent behavior 

All arguments given above may be considered as the evidence of the 
proposed statement on the equivalence of both types of the intermittent be- 
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havior. At the same time, one great difference between type-I intermittency 
with noise and eyelet intermittency taking place near the onset of the phase 
synchronization seems to exist. This difference is connected with the upper 
boundary of the intermittent behavior and this point could refute the main 
statement of this manuscript. Indeed, for type-I intermittency with noise in 
the supercritical region there is no an upper threshold (see Eq. ([1])) and the 
intermittent behavior may be (theoretically) observed for arbitrary values of 
the criticality parameter {e — Ec) > 0, although the length of the laminar 
phases may be extremely long in this case, depending on the ratio between 
the criticality parameter value and the noise intensity. Alternatively, the ex- 
istence of the boundary of the phase synchronization regime being the upper 
border of the eyelet intermittency is believed to be undeniable, since there is 
a great amount of works where the boundary of the phase synchronization 
had been observed and determined. So, this circumstance along with the 
arguments given above in Sections [TH2] involve a seeming contradiction. To 
resolve this disagreement we consider the probability P{L, e) to observe the 
turbulent phase in the time realization of the system demonstrating type-I 
intermittency with noise in the supercritical region e > Ec during the obser- 
vation interval with the length L. 

The probability to detect the turbulent phase depends on the length L of 
the observation interval and the length r of the laminar phase being realized 
at the beginning of the system behavior examination. Obviously, if r < L 
the turbulent phase is detected with the probability Pf^i = 1 and, in turn. 
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the probability for the laminar phase with length r < L to be realized is 

L 

P{t<L)= [p{T)dT, (18) 







where p(r) is given by Eq. ([5]). Otherwise, when r > L the probability to 
detect the turbulent phase is -P^^*^ = L/t, whereas the laminar phase with the 
length T takes place with the probability P{t) = p{t) dr. Correspondingly, 
the probability to observe the turbulent phase is 

L +00 

P{L,e) = I P'^^L Pir) dT+ I P^';\ p{r) dr = 



L +00 



e ne) f Le ne) (19) 
dr + I „, . dr = ^ ^ 



T{e) J T T{e 

L 



L 



L „ / L 



1 - e ne) + X r 0, 



Tie) V ' Tie) 

where r(a, z) is the incomplete gamma function, T{e) is the mean length of 
laminar phases depending on the criticality parameter e and given by Eq. ([2]) . 

The surface P{L, e) determined by the analytical expression fITI?]) and 
level curves corresponding to it are shown in Fig. |6la. It is clear that the 
probability to detect the turbulent phase for type-I intermittency with noise 
during one observation grows with the increase of the examination length L 
but decreases when the criticality parameter e — Ec is enlarged. Obviously, 
if one examines (experimentally or numerically) the system behavior in the 
time interval with the length L varying the control parameter e, one observes 
the alternation of the laminar and turbulent phases for the relatively small 

21 



1.0 2.0 3.0 
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Figure 6: The surface of probability to detect the turbulent phase in the time realization 
of the system during the observation interval with the length L and corresponding to it the 
level curves with the step AP = 0.1. The level curves = 0.1 demarcating the regions 
with the high and low probabilities to detect the turbulent phase are shown by solid lines, 
(a) The theoretical expression ([T9l) . for the simplicity the values of the control parameters 
in Eq. ^ arc taken D = 1, k = 1, Ec = 0. (b) The probability surface P{L,e) for the 
circle map with noise ((2T|) . (c) The probability surface P{L,e) for two coupled Rossler 
systems (fT3| 
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values of the e-parameter, where P{L, e) is close to one, and only the laminar 
behavior for the relatively large ones, where P{L,e) is close to zero. Having 
no information about the kind of intermittency (e.g., when the experimen- 
tal study of some system is carried out) one can suppose the presence of 
the boundary separating two different types (intermittent and steady) of dy- 
namics and, moreover, find a value 62 corresponding to the "onset" of the 
laminar behavior. Evidently, this "boundary point" would be correspond to 
the low probability P{L,e), say, e.g. Pb = P{L,e) = 0.1. In addition, one 
can perform "more careful" measurements with the increased length L' of 
the observation to determine the value of the boundary point 62 more pre- 
cisely. In this case a new value S2{L') would be obtained {P{L', e2{L')) ^ Ph), 
with it being slightly larger than the previous one. The schematic location 
of the "boundary" curve e2{L) on the plane {e,L) is shown in Fig. |6la by 
the solid line. It is clearly seen, that for the P^-level the length L grows 
extremely rapidly with the increase of the e-value. In other words, the major 
extensions of the observation interval L result in the minor corrections of 
the "boundary" point 62- Since the resources of the both experimental and 
numerical studies are always limited, some final value 62 with the maximal 
possible accuracy will be eventually found. So, despite the fact, that for 
the type-I intermittency with noise in the supercritical region the turbulent 
phases can always be observed theoretically, from the practical point of view 
(in the experimental studies or numerical calculations) the boundary point 
62 exists, above which only the laminar behavior is observed. Moreover, with 
the further development of the experimental and computational resources 
the additional studies would result only in the insufficient increase of the 
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boundary value. 

To illustrate the drawn conclusion we consider the circle map 

Xn+i = Xn + 2Q{1 — cosx„) — E + ^„,, mod 2-71 (20) 

in the interval x G [—11,71), where e is the control parameter, Q = 1.0, 
is supposed to be a delta-correlated Gaussian white noise [{^n) = 0, 
{^nC,Tn) = D6{n — m)]. If the intensity of noise D is equal to zero, the saddle- 
node bifurcation is observed in f l2U]) for e = Ec = 0, when the stable and 
unstable fixed points annihilate at a; = 0. Obviously, for the selected value 
of the control parameter Q the evolution of system (120|) in the vicinity of the 
bifurcation point may be reduced to the quadratic map 

Xn+l =xl + ^-e + ^n, (21) 

allowing an easy comparison with the results given in the previous sections. 

The intermittent behavior of type-I is observed for s < 0, whereas the 
stable fixed point takes place for e > 0. For the added stochastic force 
[D = A X 10~^) in circle map fl20l) the intermittent behavior is also observed 
in the supercritical region of the criticality parameter e, with the mean length 
T of the laminar phases and the distribution p(r) of the laminar phase lengths 
r obeying laws (jlj) and ([5]), respectively. 

The surface of the probability P{L,e) to observe the turbulent phase for 
the circle map fl20l) as well as the corresponding level curves are shown in 
Fig. To obtain this surface we have made = 2.5 x 10^ observations for 
every point taken with the steps AL = 2 x 10^, Ae = 10~^ on the parameter 
plane [e, L). The probability P{L, e) was calculated as P{L, e) = Ni{L, e)/N, 
where Ni{L,e) is the number of observations for which the turbulent phase 
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has been detected. One can see the excellent agreement between the results 
of numerical calculations and theoretical predictions (compare Fig. and 
b). 

Similarly, the analogous probability surface P{L, a) and the level curves 
shown in Fig. Elc have been calculated for two coupled Rossler systems f ll3l) 
in the vicinity of the phase synchronization boundary, where eyelet intermit- 
tency is observed. In this case = 10'^ observations have been made for 
every point to be examined, with these points being taken with the steps 
AL = 10^ and Aa = 10~^ on the plane [a, L). It is easy to see that for the 
eyelet intermittency the probability surface as well as the level curves are 
exactly the same as for type-I intermittency with noise in the supercritical 
region. As a consequence, we can draw a conclusion, that the eyelet inter- 
mittency taking place in the vicinity of the phase synchronization boundary 
and type-I intermittency with noise in the supercritical region are the same 
type of the dynamics observed under different conditions. Another conse- 
quence of the made consideration is the fact, that the phase synchronization 
boundary point can not be found absolutely exactly, since it separates the 
regions with the high and low probabilities to observe the phase slips in the 
coupled chaotic systems with the help of the experimental and computational 
resources existing at the moment of study. If someone, using a more powerful 
tools, tried to refine, say, the value of the coupling strength corresponding to 
the phase synchronization boundary reported in the earlier paper, one would 
obtain a new value being close to the previous one, but larger. Exactly the 



same situation may be found, e.g., in the work 



2l| , where two mutually cou- 



pled Rossler systems have been considered. In this work the refined boundary 
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value Eps = 0.0416 is reported with the reference to the earher work 22| . 
where the value Sps = 0.036 was given. 



5. Conclusions 

Having considered two types of the intermittent behavior, namely eyelet 
intermittency taking place in the vicinity of the phase synchronization bound- 
ary and type-I intermittency with noise, supposed hitherto to be different, 
we have shown that these effects are the same type of the dynamics observed 
under different conditions. The analytical relation between coefficients of 
the theoretical equations corresponding to both types of the intermittent 
behavior has been obtained. 

The difference between these types of the intermittent behavior is only 
in the character of the external signal. In case of the type-I intermittency 
the stochastic signal influences on the system, while in the case of eyelet 

chaotic dynamical system is used to drive the 
At the same time, the core mechanism governed 
the system behavior as well as the characteristics of the system dynamics are 
the same in both cases. 



intermittency the signal o 
response chaotic oscillatoJl 



^ Since chaotic regime has some memory in contrast to white noise, perhaps, it would 
be more appropriate to use the colored noise with a comparable memory for typc-I inter- 
mittency with noise to compare characteristics of the both types of intermittent behavior. 
At the same time, our studies show that the character of noise (such as distribution, the 
presence of memory) does not influence sufficiently on the characteristics of intermittency, 
and, therefore, the simpler model of noise may be used 
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